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1. Introduction 

In spite of the phenomenological success of the Standard Model (SM), new physics is 
expected to appear at some higher energy scale, and hopefully to provide a solution 
to the hierarchy problem, the origin of fermion masses and CP-violation, and other 
theoretical puzzles. The new physics may change the properties of the Higgs boson, 
with a substantial impact for Higgs searches in collider experiments. This is the case 
in extensions of the SM that include gauge singlet scalars, leading to Higgs boson 
decays into pairs of light neutral scalars [||] . 

Here we point out that the Higgs sector of the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) includes an axion, i.e., a pseudo-Nambu-Goldstone 
boson associated with an anomalous U{1) symmetry, which for a range of parameters 
is significantly lighter than the other scalars. The NMSSM is a well motivated 
candidate for physics beyond the Standard Model. Not only does it provide a solution 
to the hierarchy problem when combined with dynamical supersymmetry breaking. 
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but it is also free of the /i problem that plagues the minimal supersymmetric extension 
of the SM (MSSM). 

In what follows we shall investigate in some detail the Higgs spectrum and cou- 
plings in the NMSSM. The purpose of our study is twofold. First, we seek to delineate 
the region of the NMSSM parameter space consistent with the existing collider data 
which exhibits a light axion. Second, for the part of parameter space in question, 
we compute the strength of the SM-like Higgs boson coupling to axion pairs. If 
this coupling is sizable, it will have a profound effect on the collider searches for the 
Higgs boson, as the Higgs boson then decays mainly into light axion pairs, and the 
bb signature is diluted. 

The Higgs boson decay to axions persists even when the scale of supersymmetry 
is very large and the superpartners decouple. Furthermore, this phenomenon may 
occur in non-supersymmetric theories, e.g., composite Higgs models |T| or Majoron 
models [Q. 

The plan of the paper is as follows. In Section ^ we introduce our notation, 
derive the general tree-level mass matrices of the Higgs sector, and list the trilinear 
couplings between one CP-even and two CP-odd Higgs bosons in the NMSSM. In 
Section |^ we concentrate on the case of a light axion, identifying the relevant range of 
NMSSM parameters and discussing the resulting masses and mixings in the CP-odd 
scalar sector. We then derive simple analytic expressions for the Higgs spectrum 
and couplings in two cases of interest — a decoupling limit (Section ^ and large 
tan/3 (Section H). A more generic case requires a numerical study, the results being 
presented in Section p. Section ^ is reserved for our conclusions. 

2. Next-to-Minimal Supersymmetric Standard Model 

The Next-to-Minimal Supersymmetric Standard Model has the field content of the 
MSSM with the addition of a gauge-singlet chiral superfield^ S. In addition to 
the usual Yukawa-type couplings of the Higgs superfields, Hu and Hd, to the three 
generations of quark and lepton superfields, the superpotential W also includes the 
following terms involving S: 

W = XHuHdS + '^S'' . (2.1) 

We assume that the i?-parity violating terms involving quarks and leptons, and the 
dimensionful HuHa, 5*^ and S terms are forbidden by a gauge symmetry which is 
spontaneously broken above the TeV scale P]. A sector of dynamical supersymmetry 

^We shall use hatted symbols to represent chiral superfields and symbols without hats for their 
lowest (scalar) components. 
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breaking is supposed to induce masses for the squarks, sleptons and gauginos, as well 
as soft supersymmetry breaking terms involving the scalar components of Hd 
and S: 

^soft = M|Ji/„|2+M|Ji/,|2+M||5|2 + V2(mAi/>2i^d^-^5' + h.c.) . (2.2) 

The HuHd, and 5* soft terms are forbidden by the same symmetry which prevents 
the HuHd, S'^ and 5" terms in the superpotential. 

In what follows we treat the NMSSM as a low-energy effective field theory vahd 
below some scale, say in the TeV range. Therefore, we define the five mass parameters 
from Vsoit as free parameters at the electroweak scale, v ^ 246 GeV. Phenomenolog- 
ically, they are constrained by the requirement of having an electroweak asymmetric 
vacuum and a spectrum of scalars heavier than the current experimental bounds. 

The scalar potential for the Higgs sector of the NMSSM is 



XH^ia2Hd + KS''\\x''[\H, 



)\S\'' + VD + V,o{t 



with the usual D-term contributions 



Ml 
2^2 



{\Hu 



\Hd 



0%24 



(2.3) 



(2.4) 



where (Mz) is the W^-boson (Z-boson) mass. If V has a minimum where the 
vacuum expectation values (VEVs) of and Hd are aligned and non-zero, the 
electroweak symmetry is spontaneously broken. Then, in addition to the longitudinal 
W and Z {i.e. the Nambu-Goldstone modes G'^ and G°, respectively) the scalar 
spectrum includes a charged Higgs boson, H"^, and five neutral states. In general all 
five neutral states mix. However, if lm(\K.*m\m*^) <^ v^, then the scalar potential 
is approximately CP invariant. We will assume that this is the case, so that the 
mixing of the CP-cvcn neutral scalars with the CP-odd ones can be ignored. It is 
convenient to derive the spectrum using the basis where the CP-even, and 
the CP-odd, states are defined as follows: 



Hd 



^ [{v + hl- iG^) cos/? - (i/° - lAl) sin 
- G~ cos (3 + H- sin (3 



I 



Hu — 



I G+ sm 13 + H+ cos 13 \ 

\^[{'^ + K + iG"") sin (3 + {H^ + lAl) cos (3\ j 



(2.5) 
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with tan/3 = {H^) / {H^}. Notice that in order to generate masses for both up-type 
and down-type quarks one needs VEVs for both and Hd- The mxH^i(T2HdS soft 
term forces S to have a non-zero VEV, s. We choose < /3 < tt, which requires 
s > in order to minimize the first term in Eq. ( ^.31 ). 

There are several advantages of using the basis ( |2.5| ) from the beginning. First, 
notice that is rotated in the same way as v, and is exactly the linear combination of 
Hj and responsible for the masses of the W and Z gauge bosons. Consequently, 
is also the state which has trilinear couplings to W and Z pairs, and can be 
produced in association with a W or Z aX the Tevatron or LEP. In short, can 
be identified with the SM-like Higgs boson, h^. In addition, the Nambu-Goldstone 
modes and decouple from the corresponding mass matrices and need not be 
considered in our further analysis. (The basis ( p.5|) was considered also in ||^; for 
results in the more conventional basis see, e.g. [0].) 

The extremization conditions for the scalar potential allow us to replace the three 
mass-squared parameters from V^oft by the three VEVs, v sin j3, v cos j3, s: 



Mjj^ = (s^ + sin^/3) + tan/5 - ^ cos 2(3 + mxs tan (3 



2 



,^9 / 9 9 9 ^\ A/ts^ Ml „^ mxs 

Mi. = — {s' + cos^ + _ + ^ cos 2,3 + 



tan P 



\2 \^ 

M| = -— + —v^ sm2p - K^s'' + sm2(3 + m^s . (2.6) 



2s 



Therefore, the scalar masses depend on the following six unknown parameters: tan j3, 
A, K, s, mx, m^,. 

The squared-mass matrix for the CP-even scalars, H^, h^, is given by 



Ml 



( 



r + 



Ml 



V 

r cot 2(3 
A^^ — R 



—r + 



r cot 2/3 

A/t3^ + 2mxs 
v'^ sin 2(3 

-i?cot 2(3 



— R 

-Rcot 2(3 



\ 



771 



[2K^s-m^) + ^sm2(3 J 



[2.7) 



where we have defined 



(2. 



and 



R = - (Xks + nix) sin 2/? 

V 



(2.9) 
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We label the CP-even mass eigenstates in order of increasing mass, 











= U 











(2.10) 



where the 3x3 orthogonal matrix U may be obtained by diagonalizing When- 
ever there is no confusion, we shall use an alternative labelling of the Higgs boson 
mass eigenstates, in analogy to the MSSM. We shall use for the SM-like mass 
state, i.e. the one with the largest projection onto h^] for the state correspond- 
ing to the "heavy" CP-even Higgs boson of the MSSM, i.e. the one with the largest 
projection onto and H'^ for the state corresponding to the additional singlet of 
the NMSSM, i.e. the one with the largest projection onto As can be readily seen 
from ( p.7|) , in the limit of large s, the SM-like Higgs boson is identified with H^. 
However, for small values of s, h!^ can also be H2 or even ifg , depending on the other 
parameters. 

The CP-odd states, ^^,^4° have the following squared- mass matrix: 



-v{\KS-mx) (Ak + ^) sin 2/5 + 3sm«; 
The CP-odd mass eigenstates may be written in terms of the A^^ and states 



(2.11) 



A. 



cos 9 A 

— sin 6 a 



sin 9 a 

cos 6 A 




(2.12) 



with a mixing angle 9 a that satisfies 

— 4f s(Aks — m\) sin 2(3 



tan 29 a 



sm'^2p{2\KS + mx) - 2s'^{\ks + 2mx - 3m« sin 2/5) 
Finally, the charged Higgs boson has a mass 

Aks^ + 2mxs A^ 



M 



sin 2(3 



-v^ + M? 



w 



(2.13) 



(2.14) 



The vacuum defined by Eq. (|2.5| ) is indeed a viable minimum of the scalar potential 
provided all physical scalars have positive masses. Therefore, all eigenvalues of 
and M.'a have to be positive, and M^± > 0. We will analyze the constraints on 
the parameter space imposed by these conditions both numerically (Section |^) and 
analytically in certain interesting limits (Sections ^ and 

In particular, we will be concentrating on the case where one of the CP-odd 
scalars is light, and therefore the neutral CP-even scalars may decay into a pair 
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of CP-odd states. The relevant trilinear couplings in the basis ( |2.5| ) are given at 
tree-level by 



HAA 



V 

2 



y + r cot22/3 j /i° - r cot 2/5 + (^R + X^^ ) h° 



0\2 
v) 



+ 



X {k sin 2/3 + A) /i^ + Xk cos 2(3 H'^ + 2 
(^^^^K + X.hi]A.A} . 



(2.15) 



In order to compute the Higgs decay width into CP-odd scalars, one has first to 
determine the rotation matrix U and the mixing angle 6a, and then to derive the 
trilinear couplings in the mass eigenstate basis. We will perform this computation 
in the following sections. 



3. The Case of a Light Axion 

3.1 Approximate i?-symmetry 

The scalar potential V has no global continuous symmetry. However, in the limit 
where the coefficients of the trilinear terms vanish, mx,m^ 0, the potential has 
a global f/(l)_R symmetry under which the 5* charge, ys 7^ 0, is half the charge of 
HuHd- This symmetry is spontaneously broken by the VEVs of Hd and 5, so 
that apparently there is a Nambu-Goldstone boson in the spectrum. In addition, 
f/(l)ij is explicitly broken by the QCD anomaly. To see this, note that the Yukawa 
terms responsible for quark masses impose constraints on the f/(l)_R charges of the 
quarks such that the [SU{?>)cf' x U{1)r anomaly is proportional to ys- Hence the 
Nambu-Goldstone boson is in fact an axion, and there is a small contribution to its 
mass from QCD. 

Furthermore, there is another source of explicit U{1)r breaking. To see this, note 
that the form of the superpotential ( P?T| ) requires that U{1)r does not commute with 
supersymmetry, i.e. the fermion components of the Hu,Hd and S superfields have 
different f/(l)R charges than the corresponding scalars. Therefore, this U{1)r is an R- 
symmetry. Given that the gauginos carry i?-charge, it follows that ^7(1)/? is explicitly 
broken by the gaugino masses. This effect appears in the effective potential via one- 
loop diagrams with if„ ^-gaugino-Higgsino vertices. Although this contribution to the 
axion mass is larger than the contribution from the anomaly, the loop suppression 
implies that the axion is lighter than the other scalars by more than an order of 
magnitude in the limit mx, 0. 
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One has to be alert for potential confusions regarding the "axion" label used in 
this paper. This axion is not useful for solving the strong CP problem because the 
exphcit U{l)ji breaking due to gaugino masses exceeds the anomaly contribution. 
Also, the axion associated with this approximate t/(l)_R is different from the Peccei- 
Quinn axion associated with the global U{1) recovered in the k, m^^ —>■ limit. 

Another confusion may be caused by the i?-axion from the dynamical supersym- 
metry breaking sector. Typically, the models of dynamical supersymmetry breaking 
have a spontaneously broken i?-symmetry [§. The associated Nambu-Goldstone 
boson is called an i?-axion and would be massless in the absence of a source of R- 
symmetry breaking, such as a term in the superpotential required for cancelling the 
cosmological constant |p. If there is a hidden sector where supersymmetry is dy- 
namically broken, and supersymmetry breaking is mediated from this sector to the 
NMSSM via supergravity, gauge interactions, or any other mechanism, then there is 
mixing between the i?-axion and the axion discussed in this paper. However, this 
mixing is suppressed by the scale associated with supersymmetry breaking media- 
tion, and may be ignored for practical purposes. It is important however, that the 
spontaneous breaking of the /^-symmetry within the hidden sector is the source of 
gaugino masses in the NMSSM. Therefore, the existence of the i?-axion in the hidden 
sector requires a mass for the f/(l)i? axion from the NMSSM. 

3.2 Properties of the axion 

The light axion may also be identified by studying the spectrum of CP-odd states 
given in Section 2. In what follows we will expand in m\/v and vti^jv^ neglecting the 
loop effects, which is appropriate for 1 ^ rnx^^jv ^(9(10^^). The lightest CP-odd 
neutral scalar, A^, is the axion associated with the approximate f/(l)_R symmetry, 
and its mass, 

1 /2 

Ma. = VSS (m. sin^^^ + ^^^^^) + 0(</v^) , O-l) 

vanishes in the limit rrix, —>■ 0, in agreement with the arguments presented above. 
The other CP-odd neutral scalar, A2, has a mass 



Ma, ~ —JXKsm2f3 + 0{mx,m,) . (3.2) 

COSt^A 

The mixing angle < 6a < 7^/2 also has a simple form: 

tan 9a = — ^ + 0{mxjv) . (3.3) 
V sm 2p 

Notice that in the small mx, limit the axion mass Ma. depends on only four out of 
the six input parameters of the Higgs sector: s,taii i3,mx and m^. The dependence 
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on A and k drops out because they are not related to the breaking of the U{1)r 
symmetry. 

The couphngs of the axion to quarks and leptons may easily be derived by 
observing that its component has the couplings of the MSSM CP-odd scalar 
while A° does not couple to the quarks and leptons. Therefore, the couplings are 
proportional to the fermion masses within the up- and down-type sectors separately: 

— (rriu cot [3 u'y^u + uid tan (3 d'y^d) iA^ . (3.4) 

However, in contrast to the MSSM CP-odd scalar, the couplings of to down-type 
fermions are not enhanced by tan f3 because in the large tan f3 limit Eq. (|3.3|) gives 



cos^A---^. (3.5) 
s tan p 

We then see that the tan enhancement of the A^ coupling to down-type fermions 
is exactly compensated by a tan f3 suppression in the mixing angle 6a- On the other 
hand, the A^ couplings to up-type fermions are doubly suppressed by tan/5, and 
Eq. ( |3.4| ) becomes 

s (tan^ ^^^^ '^'^ d-f^d^ iA^ . (3.6) 

The phenomenological implications of this result are clear: when tan /3 3> 1 the cross- 
sections for Ai production in association with a pair of down-type fermions {e.g. bb) 
do not depend on tan /?, while the A^ branching ratios into down-type fermions are 
enhanced. 

Of special interest for Higgs phenomenology is the axion coupling to the SM-like 
Higgs boson h^, 

CkOAOAO = Iv h'A\A\ , (3.7) 

where 

ill; \ { \ — K ^\-n'}l-i\ ^ - ^ 



l + ^sin^2/5 



\[\ — K sin 2/?) + — — sin22/? + rcos 2/3 



2 



-|- Ui2 COS 2/3 I -r sin 2/3 



+ f/.3- 

V 



2k2 + ^Asin2/3( K sin 2/3 — 2/t -|- A sin 2/3 



.(3. 



Here the index i labels the mass eigenstate corresponding to /i". Although analytical 
expressions for the elements of the rotation matrix \J may be written in general 
[p!0[] , here we prefer to derive the trilinear coupling c as a series in f/s (Section ^) 
or l/tan/3 (Section ^), which should help the reader gain some insight into the 
numerical results of Section IH. 
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4. The Decoupling Limit, s ^ v 



In this Section we study the decouphng hmit of the NMSSM in which and the 
hghtest CP-even scalar, H^, are much hghter than the other scalars. This situation 
arises when the gauge-singlet VEV is large^, s ^ v, as can be seen by inspecting the 
expressions for A^^, and M^± given in Section |^. Under those circumstances, 

is simply the SM-like Higgs boson h^, and the low-energy limit of the model is 
just the SM, with the addition of one (mostly singlet) CP-odd scalar, A^. 

In order to analyze the CP-even neutral states we diagonalize by expanding 
in v/s. We keep only the leading order in itlx^k/v because we will only be interested 
in the region of parameter space where there is a light axion. The CP-even mass 
eigenstates are given by Eq. ( p.lO| ), with 

2 \ 2 / 4 \ 

Uu = l-^—{\-nsm2pf + oi^] , 



u 



12 



2k M| 



sin2/? sin4/? + C — 



U- 



13 



-sin2/3 - - 



M| cos22/3 

f2 



(A + K sin 2p) (A - 2k sin 2/3) 



A — K sin 2/3 



(SA^ - 6\k sin 2/3 + sm'^2p) - sin 2/3 cos22/3 



(4.1) 



The other elements of the orthogonal matrix U, 



21 



31 



■— sin22/?cos2/? 



s2k^ 



A^ 



+ 



-TT^ l-Tsin2/3 +0 — 



2k(A - 2/tsin2/3) Xkv"^ 

3\ 




A 



23 



22 



32 




A sin 2/? cos 2/? 
s A — 2k sin 2/3 

A sin 2/3 cos 2/3' " 
- 2k sin 2/3 , 



+ 




^Note that the dimcnsionlcss coupUngs A and k, from the superpotential of the Higgs sector are 
usuaUy assumed to be roughly of order one, i.e. we ignore cases where a hierarchy is generated 
through some kind of fine-tuning of the couphngs. 
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33 




A sin 2/3 cos 2(3" 
-2k sin 2(3 



A(A - K sin 2/?) 
2^2 



n 2" 



(4.2) 



are less important in what follows, and we only list them here for completeness. 
Notice the s/v enhancement in the last line of Eq. ( p.8[ ), which requires us to compute 
f/i3 up to an additional order in the v/s expansion. 

After computing the eigenvalues of M.\ as a power series in v/s, we find the 
following CP-even scalar masses: 



\3 

cos^2/5 ^ (A - 2fi; sin 2(3) 

2k^ 




(4.3) 



Mho 




Mh'o = V2ks + O ^ 



(4.4) 



(4.5) 



We see from Eqs. ( ^^ , ( p.ll| ) and ( |2.14| ) that to leading order in v/s, the H^, 
and scalars are degenerate, forming a weak-doublet complex scalar of mass 

^ \sm2pj 



> V 



(4.6) 



These are the familiar "heavy" Higgs bosons of the MSSM. Eq. ( [4.5|) confirms that 
the H'^ scalar is also heavy, with a mass of the order of the gauge singlet VEV s. 
Therefore, in the limit s ^ v considered in this section, the only surviving scalars 
with masses of order the electroweak scale or smaller are the SM-like Higgs boson 
and the axion A^. 

The above discussion is illustrated in Fig. where we plot the exact tree-level 
masses of the CP-even Higgs bosons (sohd hues), the CP-odd Higgs bosons (dashed 
lines) and the charged Higgs boson (dotted), as a function of s, for fixed tan/3 = 2, 
A = K = 0.5 and uix = = 1 GeV. In order to guide the eye, we have added 
shading to trace the SM-like Higgs boson h^, which from right to left is identified 
successively with H^, and H^. 

In the limit s ^ v, the condition M^o > is necessary and sufficient for the ex- 
istence of the electroweak asymmetric vacuum described by Eq. ( |2.5| ). The resulting 
constraint on the parameter space is given at tree-level by 

2k2 

instead, and the 



< -^cos^2(3+ — sin 2/5. 



(4.7) 



Alternatively, for s < i;, we see from Fig. |1| that 
requirement for positivity of the Higgs masses squared implies a lower bound of 
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Figure 1: Higgs boson spectrum as a function of s, for tan/3 = 2, \ = k = 0.5 and 
= = 1 GeV. 



s. In any way, LEP-II bounds on M^o should provide a stronger constraint on the 
parameters than the requirement of a local minimum of the potential. 

Just as an aside, notice that Eq. ( [4. 3D obeys the usual NMSSM tree-level upper 
bound on the lightest CP-even scalar [|ll|, 0, 

M|o < M| cos^2f3 + ^X^v^ sm^2p , (4.8) 




since Eq. ( [4.3|) can be equivalently rewritten as 

Mlo = Ml cos^2f5 + ^A^w^ sin22/3 - ]^X'v' ( sin 2/5 - ^ ) . (4.9) 

In this decoupling limit, hP couples to the quarks and leptons exactly like the 
Standard Model Higgs boson. Nevertheless, the decays of may be very differ- 
ent than in the SM, if Ma^ < Mh/2, since the hP — > A\A\ decay mode is then 
kinematically open. In order to assess the partial width for this decay mode 

.2„,2 / Ml 



r(/i"^A?A°) = (i_4^) , (4.10) 



we need to compute the coefficient c of the trilinear term CfiOA\A\ shown in Eq. ( |3.7| ). 
Plugging the elements of U shown in Eqs. ( |4.1| ) into Eq. ( p.8|) , after a somewhat 
tedious calculation we find a simple result: 

c= — (A-2fi:sin2/?)(A + fi:sin2/?)^ + (!:) - . (4.11) 

AKj S \ S I 
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Figure 2: Comparison between tlie analytic approximation (4.11) (dashed) and the exact 
results for |c| (solid lines), as a function of s, for tan /3 = 2, A = k = 0.5, and three different 
values of m\ = (shown, in GeV). 



This compact formula has a simple physical interpretation. In the polar coordinates 
parametrization of the axion field, the axion has only derivative couplings, suppressed 
by the axion decay constant, fA- Since the axion decay constant is given by the 
U{1)fi breaking VEV, in the s ^ f limit we find /a ~ s, while the derivatives in the 
coupling yield a factor of M^o, after taking into account the equations of motion. In 
addition, the coupling of the Higgs boson to pairs of axions has to be proportional to 
the Higgs VEV. Therefore, the amplitude for A\A\ is proportional to vM^^j . 
The amplitude has to be the same in the polar and orthogonal coordinates, which 
explains both the appearance of the Higgs boson mass M^a rather than v in the 
numerator of Eq. ([4.11|) , as well as the cancellation of the terms of order (f/s)° and 
(f/s)^. In the Appendix we provide the derivation of the Higgs coupling to axion 
pairs in the polar coordinates parametrization. 

Since A and k are expected to be of order one, the size of c is basically dictated 
by the degree of decoupling of the heavy scalars. The quality of the approximation 
is shown in Fig. |^, where we compare the prediction of the analytical formula ( [l.llj ) 
(dashed line) with the exact numerical results for the absolute value of the tree-level 
coefficient |c| (solid lines), as a function of s, for tan/5 = 2, A = k = 0.5, and three 
different values of m\ = rrif^ (shown in GeV). We see that as expected, at large s 
the analytic approximation agrees pretty well with the exact result for the smallest 
values of mx and m^. The approximation fails either at small s, or for larger m\ 
and TTift — recall that we neglected terms of order m\/v and irti^/v, while the leading 
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term in Eq. ( [4.11| ) is of order v'^/s^. The dip around s ~ 200 GeV is due to the fact 



that c changes sign, as we switch from = Hf to = Z/^- 
5. The Large tan (3 Limit 

In this section we shall obtain the coefficient c in the limit of large tan f3. To this 
end we follow the procedure from the previous Section — expand U in powers of 
1/ tan/5 and substitute the result in Eq. ( |3.8| ). However, we will not take the limit 
m\,mK 0, but instead will keep the full dependence on mA,mK, as this does not 
cause too much complication. Keeping only leading order terms, Eq. ( p.8| ) simplifies 
to 

c = -Uii + Ui2 - Ui3 2k^ - . (5.1) 

V 

The relevant mixing angles in the CP-even Higgs sector are also easy to compute — 
notice that in the large tan/5 limit, decouples and to leading order is a linear 
combination of and only. We find 

Uu = cos9h, f/i2 = 0, Ui3 = -sm9H, (5.2) 

where 

2X'^sv , , 

'""'^^=2«V-M|-m..- 
The coefficient c is then simply 

c= -A^coseH + 2/t2- sin^j^ + O ( 1 . (5.4) 

V ytan/5y 

In Fig. ^ we show the validity of the approximation (|5.4|) . In analogy to Fig. ^ 
we plot |c|, but this time versus tan/5, for a fixed value of s = 300 GeV. We see rea- 
sonable agreement between the exact result and our approximation, given the large 
cancellations between the leading order terms. We have checked that each individual 
term contributing to c is reproduced with an accuracy of a few percent, and the 
remaining difference in the total seen in the Fig. |^ is due to a large cancellation be- 
tween the first and third terms of Eq. (|5.1| ). Although the analytical approximations 
presented here are not intended to substitute the full result, they illustrate these 
leading order cancellations. It is also easy now to understand why c vanishes for a 
particular value of tan/? in Fig. |^: for sin 2/5 = X/{2k) the analytical approximation 
( PID gives c = 0. 

In conclusion of this section, we note that when in addition we take the limit of 
large s, Eq. (|5.4| ) reduces to 

^=5,0f^.^). (5.5, 



2h? tan/? 
which is also in agreement with the large tan/5 limit of Eq. ( [4.11[ ). 
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Figure 3: The same as Fig. except plotted as a function of tan /3, for fixed s = 300 
GeV. Unlike Fig. |2|, here there is a different approximation for each value of m^, as we kept 
the full dependence in this section. 

6. Numerical Results 

In this section we shall study numerically the scalar spectrum and the coupling 
( p.7|) of the Standard Model-like Higgs boson to axion pairs. There are several 
experimental (Section |6.1| ) and theoretical (Section |6.2| ) constraints on our scenario. 
Most importantly, we must account for all existing experimental bounds on the 
Higgs bosons and their superpartners. (We do not consider experimental bounds 
from other superpartner searches, since those depend on the particular framework 
of supersymmetry breaking, which we never had to specify for our Higgs sector 
analysis.) 

The Higgs structure considered here must eventually be embedded in some more 
fundamental theory, defined at a much higher scale A, possibly the Planck or the 
string scale. Hence, the low-energy Higgs sector parameters should be derived in 
terms of the parameters from the fundamental theory without excessive fine-tuning. 
It is also theoretically desirable that the couplings in the theory are free of Landau 
poles at least up to the scale A. 

Our numerical analysis in this Section is designed to address these issues. 

6.1 Collider constraints 

Let us first start with the Higgs sector. LEP is tj^ically able to rule out new particles 
with order one couplings close to its kinematic limit. For example, the charged Higgs 
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bosons can be produced in s-channel Z/7 processes. The current LEP bound is 
'mH+ ^85 GeV [|13| and we expect it to be valid in the case of the NMSSM discussed 
here as well. 

The Standard Model Higgs search at LEP, when reinterpreted as a search for 
the lightest CP-even Higgs boson of the MSSM, provides an additional constraint 
on our parameter space. The current combined LEP limit is around 113 GeV. As 
is well known however, the one-loop corrections to involving third generation 
quarks and their superpartners are positive and potentially large [Q, so that the 
LEP constraint on the tree-level Higgs boson mass Mh is much weaker. Depending on 
the stop/sbottom masses and the amount of squark mixing, the radiative corrections 



can shift the tree-level value of the Higgs mass by up to AM^ ~ 20 — 30 GeV [|15 
Hence we shall allow for tree-level Higgs boson masses as low as 80-90 GeV. Of 
course, within any given supersymmetric model framework, one can compute this 
difference exactly in terms of the parameters of the squark sector. Here we prefer to 
stay within our model-independent approach and avoid specifying a particular model 
of supersymmetry breaking and/or squark spectrum and mixing angles, as this will 
hinder the universal applicability of our results. Furthermore, the novel effect of 
h A\A\ decays is maximally operational for large masses (see, e.g. Eq.( p.5|) ), 
where the LEP bound is likely to be satisfied. 

Finally, LEP has also searched for superpartners of the SU{2) Higgs bosons, as 
part of their chargino search. Independent of their mass splitting, charged higgsinos 
lighter than ~ 72 GeV have been ruled out. For typical values of the higgsino 
mass splittings, the bound extends up to the LEP kinematic limit. The higgsino 
masses are typically of order where the n parameter, familiar from the MSSM, 
is given in our notation by 

As 

Thus the chargino mass limit constrains the product of A and s. The exact bound 
depends on the amount of mixing between the neutralinos, which again would require 
us to specify the gaugino masses within a particular model. It also depends on 
the particular framework of supersymmetry breaking, e.g. in models with low-scale 
supersymmetry breaking the bound from prompt decays of a higgsino NLSP |]17| 
from the Tevatron can be stronger. We therefore again choose to stay on a model- 
independent path, and we consider (rather conservatively) any value of ^ 100 GeV 
as allowed. 

Searches for the superpartners of the two neutral S't/(2)-singlet Higgs bosons {i.e. 
"singlinos") are quite challenging. The couplings of the singlinos to gauge bosons 
are suppressed by the neutralino mixing angles, and there are practically no limits 
coming from direct singlino production or Z-decays []TB[. The constraints on a light 
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Figure 4: Higgs boson spectrum and the absolute value of the dimensionless H^A^A^ 
coupling Cm (see text), as a function of s (in GeV) and tan/3, for A = k = 0.5 and 
nix = "T-K = 1 GeV. The horizontal (blue) dotted lines are contours of Mfjo, while the 
vertical grid of dotted red lines consists of M^o iso-mass contours. The thick (green) 
vertical line loosely marks the assumed higgsino mass bound of /i = 100 GeV, while the 
(green) shaded area contains a tachyon or is excluded from the charged Higgs search at 
LEP. The dot symbol denotes the values of s and tan/3 used in Fig. ^. The dashed line 
delineates the region where M^jo > 2M^o at tree-level. To the right of the vertically 
running dotdashed (magenta) line the SM-like CP-even Higgs boson is identified with the 
lightest CP-even mass eigenstate: = H^, while to the left of the line we have = H2 
instead. In the latter region, the corresponding coupling C211 is larger, on the order of 10% 
(see the small s region of Fig. ^). 



CP-odd scalar with small couplings to SM fermions are also very loose [^, |2^. For 
example, no significant constraints on the axion mass have been set from Z — > A^'y 
at LEP ||21|, from the Yukawa process e~^e~ —>■ ffA^ ||2^, from direct production 
through gluon fusion at the Tevatron |]23|, from fits to the electroweak data [^], 
or from meson decays [||]. The LEP searches for Z* h^A^ A^A^A^ within 
the MSSM can be reinterpreted as axion searches, but the limits are diluted 
because the axion coupling to the Zh^ in our case is smaller by a factor of cos 6 a 
[see Eq. (p.5|)]. The relevant lower bounds on Ma currently come from beam dump 



experiments ||2^, in the MeV range, and from star cooling rates. Ma ^0.2 MeV [p^j . 

Having summarized the relevant collider constraints, we now present our exact 
numerical results for c. In Fig. ^ we show contours of the masses of the SM-like Higgs 
boson, M/jO, [horizontal dotted (blue) lines] and the lightest CP-odd Higgs boson. 
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My^o, [vertical dotted (red) lines] as a function of s and tan/?, for fixed \ = k = 0.5 
and m\ = = 1 GeV. The thick (green) vertical line denotes the assumed higgsino 
mass bound of = 100 GeV and the region to its left is disfavored. Inside the (green) 
shaded area there is either a tachyon in the spectrum, or the charged Higgs mass is 
below the experimental bound. While it is in principle possible that the radiative 
corrections to the Higgs boson masses, neglected here, may shrink, or even eliminate 
this excluded region, we do not find this latter part of parameter space particularly 
attractive, since it is associated with very light higgsinos. 

Throughout most of the parameter space shown in Fig. ^, the SM-like Higgs 
boson is the lightest CP-even mass eigenstate: = Hi. However, at small values of 
s, and H'^ become light as well (see Fig. |l|). The dotdashed (magenta) line going 
vertically marks the point at which the SM-like Higgs boson changes its identity — 
from /i° = (to the right of the line) to h° = (to the left). 

In Fig. 1^ we only show the masses of the and Higgs bosons, since we 
are interested in the phenomenology of the A^A^ decay. The first important 

result seen in Fig. ^ is the location of the region of parameter space, where this decay 
is open. The dashed line in Fig. § delineates the relevant part of parameter space 
where Mj^o > 2M^o at tree-level. The location of this region is easily understood. 
Recall that M40 scales with -y/i, as evident from the Figure, as well as Eq. (|3.1| ). 
Then notice that unlike the case of the MSSM, here the lightest CP-even Higgs 
boson mass decreases with tan /?, so that tan /? values as low as 1 — 3 are possible. 
Combining these two observations, we easily see that the /i° — >• A^A^ decay is most 
likely to be open at small tan/3 and small s. 

Our second main result shown in Fig. ^ is the strength of the Higgs to axion 
coupling. We show contours of the absolute value of the dimensionless coupling 
I Cm I, which we define in analogy to Eq. (p.7|) : 

^HOAOAl = I E c.^'c H^A'^^Al . (6.2) 

As we already explained above, to the right of the vertical dotdashed magenta line. 
Cm is identical to the coefficient c defined in Eq. ( p.7|) . We see that in the region of 
parameter space, free of any experimental constraints, c can be as large as 0.05. We 
also see the possibility of exact cancellation and vanishing c. Indeed, in the limit of 
large s, our leading order approximation Eq. ( [4.11| ) vanishes for tan jS = 2 + ^/3 ~ 3.7, 
in reasonable agreement with Fig. ^. In the small s region, where h'^ = H2, the 
corresponding C211 can be even larger, on the order of 10% or more. 

In summary, our main conclusion from Fig. ^ is that for fixed A and k, values of 
c are typically maximized at small s and small tan f5 — exactly in the spot where the 
— > A^A^ channel is most likely to be open. This could also have been anticipated 
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Figure 5: The same as Fig. ^, but as a function of A and k, for s = 300 GeV, tan (3 = 2 
and mx = m-^ = 1 GeV. The dot symbol shows the values of A and k used in Fig. ^ 

from the approximate scaling c ~ M^o, which we found in the two limiting cases 
considered earlier. 

We now turn our attention to the dependence on the other two main parameters: 
A and k. To this end, we fix s = 300 GeV and tan /? = 2, and show the Higgs spectrum 
and |ciii| in Fig. ^ as a function of A and k, again with fixed mx = rrii^ = 1 GeV. 
The (green) shaded area is again excluded because of tachyonic or light charged 
Higgs states, and the (green) vertical line denotes the assumed higgsino bound /i = 
100 GeV, with the region to the left of it being disfavored. This time the SM-like 
Higgs boson is unambiguously identified as = H^. On the other hand, the axion 
mass is completely fixed in terms of mx, m^, s and tan/5 [recall the discussion after 
Eq. ( p.3|) ]. For the values of the parameters considered here, M^o = 33.72 GeV 
and its fractional variation throughout the whole Figure is less than 1 part in 10^. 
The dashed line, depicting the region M^o > 2M40, is therefore coincident with the 
contour of M^o = 67.44 GeV. 

We see from Fig. ^ that large ratios of A/k are disfavored. This can be easily 
understood from Eq. (^4.3| ) — the second term gives a large negative contribution 
to the lightest CP-even Higgs boson mass, which results in M^o < 0. Nevertheless, 
there is a significant allowed region with naturally large couplings A and k. Now 
the lightest CP-even Higgs boson mass can be larger than the LEP limit already at 
tree-level. The coefficient c is again maximized in the region with the largest M^o, 
and again values of c ~ 0.05 are possible. 



18 





Given that the small axion mass allows for the — > A^A^l decays in principle, 
and that the coefficient c can be sizable (see Figs. | and it is interesting to quantify 
how the branching fractions of the SM-like Higgs boson, h^, are affected. To this 
end, one must use the physical mass Mjp^^, as the partial widths for A^A^ and 
h'^ — i> bb have the opposite dependence on M^^^ . So, by using the (smaller) tree- level 
mass, one would be overestimating T{h^ ^i^i) ^ind underestimating T{h^ bb). 
This is why in Fig. ^ we show the branching ratio for A^A^ decays 

as a function of the physical Higgs mass M^^^ and the value of the coefficient c, 
assuming M^o -C M^^^ and M^o ^ Mf^^^ . Here we have taken Fsm as the width 
of the SM Higgs boson. Fig. |^ should be interpreted as follows. After fixing the 
fundamental input parameters A, k, m^, m^, s and tan/?, one should compute the 
physical Higgs boson mass Mf^^^ in terms of the parameters of the squark sector 
derived within a given model. Then the coefficient c can be read off from Figs. ^ and 
^, and the corresponding branching fraction B{h^ — > A\A^ for the resulting values 
of c and M^o^^ is given in Fig. |^. 

From Fig. ^ we can see the potential importance of the A^A^ decay mode. 
In the theoretically preferred mass range Mlo""" ^ 160 - 180 GeV of the NMSSM, 
this decay mode completely dominates for c ^0.1. We have seen that such values 
of c are in principle possible, although in limited regions of parameter space. For 
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the more typical values of c ~ 0.02 — 0.05, the Higgs boson branching fraction into 
axions is comparable, if not larger, than B{h — >■ bb), and may still significantly affect 
the Higgs boson collider searches utilizing the bb mode. 

6.2 Theoretical prejudice 

The theoretical constraints discussed in this subsection are less robust, since one must 
define a framework. For example, when requiring the absence of Landau poles for 
the Higgs couplings A and n, one must specify the high-energy scale A, at which the 
effective theory description of the NMSSM breaks down. Naturally, the constraints 
would be stronger, if A is identified with the GUT or Planck scale, rather than some 
intermediate scale, related to the (mediation of) supersymmetry breaking. 

The Renormalization Group equations (RGB's) for the Yukawa couplings of the 
NMSSM have been extensively studied in relation to the question of the absolute 
upper limit on the lightest CP-even Higgs mass p9|, EO]. If A is identified with 



the GUT scale, the perturbativity requirement implies A ^0.87 and k ^0.63 W2 



If, however, A is identified with an intermediate scale, as in gauge- 



mediated models [TO, then the bounds can be somewhat relaxed, e.g. k ^ 1.36 



for A ~ 50 — 100 TeV [|T^. We see that these constraints still leave a lot of the 
available parameter space in Figs. ^ and |^, where the new decay A^A^ is 

possible, and its branching ratio is sizable. 

One may also wonder if the low values of and needed in order to suppress 
the axion mass Mai (see Eq. ( |3.1|) ) can appear naturally, without any significant fine- 
tuning. So far we have adopted a low-energy point of view and never specified the 
particular model framework for supersymmetry breaking and its communication to 
the NMSSM sector. The size of nix and will depend on two factors: the boundary 
conditions at the high-energy scale A, where the soft supersymmetry breaking pa- 
rameters are generated, and second, the amount of (logarithmic) RGE running from 
that scale down to the electroweak scale v. Notice that due to the singlet nature 
of 5", the one-loop beta function for only depends on and m^. In the ?7(1)r 
limit both of these two parameters start out small at the high-energy scale A, and 
the generated value for at the weak scale is typically also rather small. On the 
other hand, the one-loop beta function for depends on the gaugino masses Mi 
and M2, as well as the rest of the A-term parameters, which can be relatively large 
already at the scale A. Thus the induced value for at the weak scale can be much 
larger than m^. However, a closer inspection of Eq. ( |3.1| ) reveals that the axion mass 
is usually more sensitive to ttIk rather than nix. For example, in the large tan/3 limit 
we have cos 6^^ ~ 1/ tan (3, hence 

Mlo ~ O {ni^s) + O j ~ O M . 



20 



Without a particular model, it is difficult to be more specific at this point. We shall 
therefore leave this question open for future studies. Let us only point out that there 
exist at least two very well motivated frameworks, in which the boundary conditions 
for mx and are zero at the scale A — gauge mediation ^] and gaugino 
mediation ||35 |. 



7. Conclusions 

We have discussed a limit of the NMSSM where a light CP-odd scalar (axion) is 
present in the Higgs spectrum. The light axion appears as a result of an approxi- 
mate global U{1)r symmetry of the scalar potential, which is spontaneously broken. 
We found that the mass of the axion is proportional to the soft breaking trilinear 
couplings mx and m«;, and if those are in the GeV range, the axion can easily be 
lighter than half the SM-like Higgs boson mass, M^. In those cases, we computed 
the h^A^A^ coupling and found that it can have a direct impact on phenomenology, 
as it can substantially modify the SM-like Higgs boson collider signatures. 

In conclusion, we feel that our results fill a major gap in the extensive liter- 
ature on the Higgs sector of the NMSSM. Previous NMSSM studies have mostly 
concentrated on setting absolute upper limits on the SM-like Higgs boson mass PB | 



(which is of primary interest for the production of in collider experiments) or 
the related singlino phenomenology |jT8|, However, the case of a light CP-odd 
axion considered here has largely been overlooked. In light of the interesting phe- 
nomenological implications of the scenario presented here, and the symmetry reasons 
behind its motivation, it is worth pursuing the case of a light CP-odd scalar in future 
phenomenological studies. 
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Appendix: Higgs Coupling to Axions in Polar Field Coordi- 
nates 

In this Appendix we derive the Higgs boson coupling to axion pairs using the polar 
field coordinates. Although the physical results are independent of the parametriza- 
tion of the degrees of freedom, and the orthogonal field coordinates used in section 
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2 [see Eq. ( p.5| )] are convenient enough, the polar field coordinates are particularly 
appropriate for describing the axion in the limit where we neglect its mass. 
Consider the following parametrization: 



f ^[{v + hi) cos/3 - sin/3] e^(^St^'^/5-^?°)^ ' 



V 



- G- cos 13 + R- sin (3 



H„ 



S 



f G+ sin (3 + H+ cos f3 

^j^[{v + hi) sin/3 + cos/3] e*(^S cot/3+G«)/^ 



s + h\ 



(A.i; 



With this parametrization, it is straightforward to derive the masses and mixing 
angle of the CP-odd states obtained in section 3.2. 

More importantly, the scalar potential Eq. ( p. 51) does not involve the axion field 
in the mx,K. limit. If we neglect the axion mass, then only the kinetic terms give 
rise to axion couplings. This is the usual statement that a Nambu-Goldstone boson 
has derivative couplings in the polar field coordinate parametrization. The kinetic 
terms for Hd and S include canonically normalized kinetic terms for the physical 
states as well as derivative couplings of CP-odd scalar pairs with CP-even scalars. 
The SM-like Higgs boson coupling to axion pairs can be easily derived: 



c'v 
2M! 



(A.2) 



where the dimensionless parameter c' is given by 
M? ^ 



V 



{Un + 2Uu cot 2/3) cos'Oa + -f/is 

sin 6 a 

s 



(A.3) 



In the s ^ V limit, the elements of the matrix U that rotates the CP-even scalars to 
the mass eigenbasis have been computed in section 4. Using Eq. ( [4.1| ) we find a very 
simple result: 



c' = -2c[l + 0{v/s,mxjs)] 



(A.4) 



where c is given by Eq. ( [4.11| ). The last step of this computation is to check that 
the decay width for h^ —>■ A^A^ is the same as the one obtained using orthogonal 
coordinates [see Eq. ( |4.1CI| )]. This straigthforward exercise provides the explanation 
for the M| / dependence of c in Eq. ( |4.11| ) . 
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